
  #include "math.h"
  void powr(x,y,n,u,v)
  int n;
  double x,y,*u,*v;
  { double r,q;
    q=atan2(y,x);
    r=sqrt(x*x+y*y);
    if (r+1.0!=1.0)
      { r=n*log(r); r=exp(r);}
    *u=r*cos(n*q); *v=r*sin(n*q);
    return;
  }

